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We have developed a rigorous quantum model of spontaneous parametric down- conversion in a 
nonlinear ID photonic-band-gap structure based upon expansion of the field into monochromatic 
plane waves. The model provides a two-photon amplitude of a created photon pair. The spectra of 
the signal and idler fields, their intensity profiles in the time domain, as well as the coincidence- count 
interference pattern in a Hong-Ou-Mandel interferometer are determined both for cw and pulsed 
pumping regimes in terms of the two-photon amplitude. A broad range of parameters characterizing 
the emitted down-converted fields can be used. As an example, a structure composed of 49 layers 
of GaN/AlN is analyzed as a suitable source of photon pairs having high efficiency. 

PACS numbers: 42.50.Dv 



I. INTRODUCTION 



More than 20 years ago. Hong, Ou, and Mandel showed 
experimentally that mutually strongly quantum corre- 
lated (entangled) photon pairs can be emitted in the 
nonlinear process of parametric down-conversion [l|,l2| at 
the single-photon level. For the occurrence of entangled 
photon pairs, the spontaneous character of the process is 
important. Entanglement of two photons comprising a 
photon pair might occur for various physical quantities 
like frequencies, emission angles (wave- vectors), or polar- 
izations. Perhaps most interestingly, entangled photon 
pairs manifest themselves in time domain, where both 
photons are detected within a relatively narrow time win- 
dow. This is a direct consequence of the 'point' character 
of the emission of a photon pair in the time domain. The 
width of the time window characterized by an entangle- 
ment time is typically on the order of several hundreds 
of fs, and is experimentally observable using Hong-Ou- 
Mandel interferometer. 

In the time that has intervened since the original pre- 
dictions, entangled photon pairs have been used for nu- 
merous experiments demonstrating both fundamental as- 
pects of their physical properties 3] and their potential 
for applications. These also include tests of Bell inequal- 
ities PI , quantum teleportation H, the generation of 
Greenberger-Horne-Zeilinger states |3 or quantum com- 
putation [Q] that can exploit photon pairs. Quantum 
cryptography with photon pairs Q niay be considered 
the most important application. We mention that appli- 
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cations to metrology have also been suggested Q. 

Most researchers in the field have focused their atten- 
tion for the most part on bulk nonlinear crystals, pumped 
by intense laser beams in type I and type H configu- 
rations. Bright sources of polarization-entangled pho- 
ton pairs have been fabricated using two nonlinear typ e 
I crystals, mutually rotated by 90 degrees ^ [ifl Um- 
The nonlinear processes in an optical cavity have also 
been used to enhance the photon-pair generation rate 
|l2l |. In addition, periodically-poled materials may be 
used to increase the photon-pair generation rate in ma- 
terials where phase matching is not naturally available 
[l3|. At present, even sources of photon pairs pumped 
by laser diodes have been developed [iJI- 

On the fabrication front, techniques of structures com- 
posed of nonlinear thin layers (of width of several tens or 
hundreds of nm) have developed to the point that useful 
structures wit h p re-defined properties may be achieved 
rather easi ly IJ SI ■ These nonlinear photonic-band-gap 
structures lla llTl Il8l Il9l | are very promising as sources 
of photon pairs, as has recently been shown in |2(1 l2ll|. 
Despite the small amount of nonlinear material embed- 
ded inside them, they can generate photon pairs with 
relatively high efficiency thanks to the constructive in- 
terference that involves both the pump and the down- 
converted fields, their spatial inhomogeneity nothwith- 
standing. The enhancement of the photon-pair genera- 
tion rate has been predicted to be several hundreds and 
even thousands of times larger than photon-pair genera- 
tion rates in nonlinear, bulk materials 21]. Moreover, the 
spectral and spatial characteristics of the down-converted 
fields depend on details of the structure, a fact that might 
be used to control the process, at least to some extent. 
For example, photon pairs with very narrow spectra may 
be obtained from suitable structures. We note that along 



the same vein, four-wave mixing in photonic-band-gap 
nonhnear fibers is also promising as a modern source of 
photon pairs |2ll2^. 

In this paper, we present a quantum model of photon- 
pair emission in a nonlinear, one-dimensional photonic- 
band-gap structure, based upon a perturbative solution 
of the Schrodinger equation. This model extends those 
developed for bulk nonlinear materials in ^24 .25.- ■26. . .271 . 
It has been shown in J23| that this approach is compat- 
ible with models based on methods of classical nonlin- 
ear optics (multiple-scale spatial and temporal expansion 
methods) with a specific kind of stochastic averaging con- 
cerning emitted spectra. The quantum model, however, 
provides a complete description of the nonlinear process. 

The paper is organized as follows. In Sec. II, the 
model is presented in three steps. Description of a quan- 
tum field in a layered medium in Subscc. II A is followed 
by the description of nonlinear quantum interactions in 
Subsec. IIB. In Subsec. IIC measurable characteristics 
of the generated down-converted fields are determined. 
The behavior of physical quantities characterizing a pho- 
ton pair are discussed in Sec. Ill, both for cw and pulsed 
pumping regimes. Sec. IV contains our conclusions. 
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FIG. 1: Scheme of the considered nonlinear one-dimensional 
layered structure: nl , n\ , and Up denote indices of refrac- 
tion of signal, idler, and pump fields in Ith layer with the 
length Li, d''' stands for a nonlinear coefficient of Ith layer 
and zi are positions of the boundaries. 
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II. DESCRIPTION OF SPONTANEOUS 

PARAMETRIC DOWN-CONVERSION IN A 

ONE-DIMENSIONAL, NONLINEAR 

PHOTONIC-BAND-GAP STRUCTURE 



We consider a stack of A^ nonlinear layers in which 
spontaneous parametric down-conversion may occur. A 
sketch of the system is shown in Fig. ^ The Ith layer 
begins at z = Z(_i and ends at z = zj, its length is de- 
noted as Li, I = 1, ... ,N. Linear indices of refraction of 
pump, signal, and idler fields in the Ith. layer are denoted 
as rip , ni' , and nl , respectively. Symbols ni ' , nl , 

and Up {lis J "-i 1 and rip ) mean indices of 
refraction in front of (beyond) the sample. The symbol 
d*-'' is used for the nonlinear tensor of Ith layer; symbols 



,(0 1.(1) 



and fcj- '' stand for pump-, signal-, and idler- 



field wave-vectors in the Ith layer. 



Optical fields in a one-dimensional 
photonic-band-gap structure 



The structure is pumped by a classical strong pump 
field that propagates under the angle "dp with respect to 
z axis (see Fig.|2l. Its wave- vector kp lyes in the yz plane, 
whereas its electric-field amplitude Ep^^ is perpendicular 
to the wave- vector kp. The pump-field positive- frequency 
electric-field amplitude at the input of the structure at 
z — zq IS denoted as Ep^, {zo,t) and can be conveniently 
described using a positive-frequency electric-field ampli- 



FIG. 2: Scheme of the used coordinate system for the descrip- 
tion of pump, signal, and idler fields. A field with wave- vector 
k (lying in yz plane) and polarization vector e (perpendicular 
to k) propagates at an angle i9 with respect to z axis. The 
field is linearly polarized and the polarization vector e makes 
an angle ip with respect to x axis. The polarization vector 
exE of TE wave is parallel to x axis whereas the polarization 
vector eTM of TM wave lyes in yz plane. 

tude spectrum Ep^ (zo,^p) determined as follows: 

eI,+\zo,^p) ^ ^= diEW(zo,t)exp(ic^pi). (I) 

If we assume a Gaussian time profile and a linear po- 
larization of the incident electric-field amplitude Ep^ in 
the direction rotated by an angle tpp with respect to an 
incident TE-wave polarization direction ep^xE i^pF,TM 
denotes an incident TM-wave polarization direction, see 
Fig. [2Jl we have: 
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where ^p is the pump-pulse amplitude, Tp pulse duration, 
LUp central frequency, and Up denotes a chirp parameter 

of the pulse. The spectrum Ep^, {zo,ujp) determined by 
Eq. (Q) is given as: 



E(+)( 
PF y 



Zq, LUp 



[epj,,TE cos{(pp) + epp^TM sin{ipp)] 



xCp 



^2(1 + iflp) 



exp 



4(1 + mp) 



K-^^)' 



= e(,+|j,^(zo, ^p) + E^+^_j,^,,(zo, ^p). (3) 



On the other hand, the foUowing spectrum corresponds 
to cw pumping: 



E],+>{zo,u;p) = [ 



GpF.TE COs((pp) + Bp^^TM Sin((/Jp)] 

xep<5K-c.°); (4) 



5 means a Dirac delta function. 

The pump field incident on the structure is scattered 
at each boundary inside the structure, so as to achieve 
a certain profile along the z axis, provided that the non- 
linear interaction does not lead to pump-field depletion. 
Scattering of the pump field is conveniently described 
using its decomposition into monochromatic waves. The 

positive-frequency electric-field amplitude Ep^a{z,ujp) of 
a monochromatic component at frequency ujp with polar- 
ization a {a — TE, TM) can be written as follows |28i |: 



E(,|j(z,Wp) = rect_oo,zo(^) 

X 4"J, Ja;p)eW^ J^p) exp(zk(«) (z - zo)) 

+4«J^Jc.p)eW_Ja;p) exp(~*k(«) (z - zo))' 
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X [Al!ljLUp)el!ljcOp) exp(*k« (z - z,_i)) 
+41aK)4i«K) cxp(-^kW (^ - zi^i)) 

X A(^+i) (u;p)e(^X'^ (-.) exp(*k(^/^) (. - ^a^)) 
+4^+i)(a;p)e(^ii)K) exp(-*kf/i)(z - z^.)) 
a = T^, TM; 



(5) 



the function rectz^^2^(z) equals one for Za < z < Zb and 
is zero otherwise. Polarization vectors of a-waves in Zth 
layer are denoted as Gpl^a and ep^^Q for forward- and 
backward-propagating fields with respect to z axis, re- 
spectively, and they are frequency dependent. The sym- 
bol kp/z denotes a z component of the pump-field wave- 
vector kp in Ith layer and is determined by the expres- 



k« =kWcos(,9W), 



(6) 
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where the wave- vector kp propagates in the Ith layer 
under the angle -(^p with respect to z axis. The angles 
•dp fulfill the Snell law at the boundaries, i.e. 

4') sin(4')) = 4'+i) sin(i9(,'+i)), l = 0,...,N; (7) 
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The symbol ApJ,a{iLJp) is identified as an electric-field 
amplitude of the pump field at frequency ujp with po- 
larization a incident on the structure from the left- 
hand side, whereas the symbol Ap^^a (wp) describes an 
electric-field amplitude of the pump field at frequency 
ujp with polarization a incident from the right-hand side. 

The remaining amplitudes Apl^a{ujp) and Apl^a{u}p) are 
determined using relations at the boundaries and free- 
field propagation inside the layers: 



i(i) 



,(0) 



App,ai^p) \ _ ^(0)/, , \ ( App^aiuJp) \ 



PF,a{l^p) \ _ ^(l)( y-pil)),. \ ( App,a{uJp) \ 



PB,a y^p 



i^p) 



_ Apg^a{l^p) _ 

a = TE,TM; l = l,...,N. (8) 



We note that the coefficients Apl^a and Ap^^a for I = 
1, . . . , TV describe the corresponding electric-field ampli- 
tudes at the beginning of /th layer. 

Assuming TE and TM waves, the boundary transfer 
matrices T^j^^ and T^tm have the form: 
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2 l^ 1 - /(') (^p)g« K) 1 + fP (u:p)c/i> (cUp) ^ 
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1 ( fp {^p) + 9p i^p) fp i^p) - gp i^p) 

2 V fp^ i^p) - ai'^ i^p) fp^ H) + ai'^ H) 



; = o, 



.N: 



(9) 



fj>^ = cos(4'^)/cos(4'^'^) and g^P = n'p^ /n'p+^l The 
free-field propagation matrices Vp (cLip) can be written 



as: 







' ' V exp(-*4',U0 

l^l,...,N. 



(10) 



The positive-frequency electric-field operators 
^^'{z^t) and E^{z,t) for the signal and idler 
fields can be decomposed into TE- and TM-wave 
contributions Ey'^^{z,t), Ei^"^ j^^{z,t), E]'tji^{z,t), and 
E^ j?^(z,t) and expressed as follows l29j : 
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Em.T_E(^''^m) +Em.TM(^'^m) 



m — s,i. 



(11) 



The permittivity of vacuum is denoted as eg; c means 
speed of light in vacuum, h is the reduced Planck con- 
stant, and B denotes the area of the transverse profile 
of a beam. The symbols e,„p_Q,(z,cjm) and emB,a{z,ujm) 
mean polarization vectors of mode m of a wave propagat- 
ing forward and backward with respect to the z axis. An- 
nihilation operators of signal [asp,aiz,uJs), asB.a{z,ujs)] 
and idler [0,^^0,(2,^^), aig,a{z,uji)] photons in a wave in- 
troduced in Eq. (|ll(l can be expressed in the same way 
as the pump-field amplitude Ep,a{z,ujp): 

— reCt_oQ 2:0 v-^J 

+ amLa('^™)4°L.a(^m)exp(-ik(^)^(z - Zq)) 
N 

-H^rect^,_,,^,(^) 
1=1 

X [a'2lA'^m)e2l,aiuJr,i) exp(ik(^^^(2 - zi^i)) 

+ a^B.a('^™)emB,a(^m) exp(-ik(^^,^ (^ - Zl-l)) 

-|-rectz„,oo(^) 

X '&itta^ Kn)e(f,+^) (c.„) exp(zk(,^+i) (z - z^)) 

m = s,i- a = TE,TM. (12) 
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The polarization vectors einp.a and emg.a give the po- 
larization directions of mode m with a wave in rth layer 
propagating forward and backward, respectively, whereas 
km,2 is a z component of the wave vector km of this 
mode: 



kW =k«cos(^«), m = s,z. 



(13) 



,(0 



The angle dm characterizes the direction of propagation 
of mode m with respect to the z axis. The angles -dm are 
given by the Snell law at the boundaries, i.e. 

nWsin(^«) = re^^Bin(^^i)), 

TO = s,i; / = 0, ...,iV; (14) 



"dm — dm, where dm stands for the angle of incidence of 
mode TO. 

The operators amp,a{i^m) and a„ii3,a(wm) obey the fol- 
lowing commutation relations: 

[ai'^pJ^m),alll\,{uj'm)] = Sa,a'Sm.,m'Sl,l'S{uJm-Uj'm): 



^(Um), a I ,(llJ, 



m'.a' \ mJ 



4)] = 0, 

-'m)] — Sa,a'^m,m'Sl,l'S{uJm — ^m), 

4)] = 0, 

4)] = 0: 

4)] = 0; 

m,m'^s,i; a^TE,TM; 1,1' ^ 0, . . . , N + I. (15) 

The photonic-band-gap structure imposes the fol- 
lowing relations among the operators amp a{^m) and 
a(0 



[amB,a(^™)'am^,a'('^" 

[a?Jp,o,{uJm),allJ^^^,i^'„ 

[al^pA^"^)'^m'[a'(^'m, 



.q(cl;„i) acting in Ith layer: 






a{^m) 



r(") (a 



a(") 



a\(^m) 



(0) 

a(0 



^77), R .a v'^m/ 



+ 1)/,, N I - ^m,a' rr, \ (I) , . I , 



m^s,i] a^TE,TM; l^l,...,N. (16) 

Transfer matrices Tmja at the boundaries and free-field 
propagation matrices Vm for m = s,i and a = T_E, TM 
are defined in the same way as those given in Eqs. 
and pU|l for the pump-field amplitudes. 



B. Nonlinear interaction inside the 
photonic-band-gap structure 

The Hamiltonian H{t) describing spontaneous para- 
metric down-conversion in a nonlinear medium of volume 
V at time t can be written as: 

H{t) - eo / dr 
Jv 

d(r) : [E(+)(r,i)E(-)(r,i)Ei-)(r,i)+h.c.] , (17) 

where d denotes a third-order tensor of nonlinear coef- 
ficients, the symbol : is shorthand of the tensor d with 
respect to its three indices, and h.c. stands for a hermi- 
tian conjugated term. The negative-frequency electric- 
field operators Em for m = s,i have been introduced 
in Eq. I^T)) (Em ' = Em ). Decomposing the interact- 
ing fields into TE and TM waves, and using the inverse 
Fourier transformation in Eq. p7|l . we arrive at: 



H{t) = ^ / dz r dujp 

V27r Jo Jo JQ 



diOf, I duJi 
'0 



E 

a,p,i=TE,TM 



d{z) : EW(z,c.^)Elj(z,w.)Ey(z,^0 



(18) 



Integrations over the variables a; and y in Eq. p7|l im- 
pose conditions for x and y component of wave- vectors; 



S{]s.s^x + K,x) [^p,x = is assumed] and S^ks^y + k^^j, - 
kp y). The latter (5-function provides the foUowing rela- 
tion between the angles -di ' and -dl ' of modes s and i in 
the ^th layer: 
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(0 



aresm 



^sin(^«)-i^sin(^«) 

.,N+1. 



(19) 



Solution of the Schrodinger equation to first order in 
nonlinear perturbation together with the assumption of 
incident vacuum state |vac) in signal and idler fields pro- 
vides the output state |''/')s"/3 i 7 °^ signal field with j3 po- 
larization and idler field with 7 polarization in the form: 
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s,/3,z,7 



I vac) 



2\/2ttc 



N 
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E 
1=1 



dujp I duJs^/uTs 
"'0 



X / duj^^^Y. Y. Y. Y. 

f* 1=1 m=PF,PB 11 = Sp,SB 0=ip,iB 

X Y. dW : e« Jc.,)ei'yc.,)eW (^,) 

a=TE,TM 
X^m,aK)iiexp 



+ ^(i^W_Al')_i,(0)i^ 



l(X«-XW-if«)L, 



xai';j(c..)a»t(^,)|vac). 



5{ujp - LUs - idi) 

(20) 



The wave- vectors X*^'^ introduced in Eq. H2U|I are defined 



as K 



(0 



k^.'l and K^^ 



3F J. -2 J-E 

(0 



= — k^ ^ for j = p, s, i. 

and drriB.a for a waves in mode 



The operators amp, a 
m in Ith layer can be expressed in terms of the operators 

amp. a and amg.a- These relations can be, e.g., written 
in the form: 



«(') (,, \ 






^/{Sni,a)ll{^m) ''{Sm^a)l2i^m) / {Sm^a)ll{i^m) 







1 



m = s,i; a — TE, TM\ I — 1, . 



,N. 



(21) 



The matrix Sm,a used in Eq. 121|) describes the propa- 
gation of a wave in field m through the whole structure: 



m — s,i; a = TE, TM. 



(22) 



Similarly, the pump-field amplitudes Apl^a and Apl^a 
for a waves in Ith layer can be determined from the ampli- 



tudes of the incident fields Apl^a and Ap^^a as follows: 



App,a(iOp) \ _ 
ApB,a\'^p) I 



Yl[v'^'\u^p)T}A'\u^p) 
3=1 







\ —{SpA2l{^p)/{SpA22{^p) 1/(5^^0)22 (Wp) 

y ApB,a y^p) J 

a^TE,TM; l^l,...,N. (23) 

The matrix Sp^a describes the propagation of a classical 
pump field with polarization a through the whole struc- 
ture, i.e. 

1 

SpA^p) = Tp^HKu^p) n [v^p'\udp)T;;A'\u:p)\ . (24) 

We note that the expression in Eq. H2(J|I for the output 
state \4')°s^p i ^ including relations written in Eqs. (|21I24() 
can be formally recast into a compact form using the so- 
called left-to-right ($(+') and right-to-left ($(^)) modes 
introduced in the classical electro-magnetic theory of lay- 
ered structures (for details, see |21|). Fields exiting the 
structure at z = z^v are described by $(+^(z) functions, 
whereas $'•^-'(2) functions are appropriate for fields exit- 
ing at z = zq. In classical theory, this corresponds to the 
picture in which every nonlinear layer is a source (emit- 
ting 'dipole') of photon pairs P0| and the fields leaving 
the structure are given as a sum of contributions from all 
layers. 

We assume that the outgoing signal (idler) field is de- 
tected using an analyzer with the polarization that forms 
an angle ips {(pi) with respect to the TE-wave polariza- 
tion direction (see Fig. EJ . In order to get the right wave- 
function, we transform the operators of the outgoing sig- 
nal (idler) fields into the basis with the polarization vec- 
tors given by angles tps {(fii) and ips+Tr/2 (}pi + Tr/2) using 
the following formulas: 

/ "■m.F^TE \ _ 

{N+l) 



a 



'rap 



,TM . 




[N+l) 
(N+l) 



(0) 

•mB.^m 

d^°^ 

"lB,'P,n+Tr/2 , 



(25) 



Substituting relations in Eqs. ||2U, if^ . and if^ into 
the expression \ip)°s J^E,i,TE + W°sJ^E,^,TM+\i')°sJ'M,^,TE+ 
IV')s"tm i TM determined using Eq. (^01) . we arrive at the 



output state |i/')™^^,i,<^i describing a signal photon po- 
larized along the angle t/jj and an idler photon polarized 
along the angle ipi : 



\ out 



vac) 



dw, / duj. 






*■ ■" a) ,„. vac) 






.(26) 



The function (^""(w^jWi) introduced in Eq. l|^ has the 
meaning of probability amplitude of having a signal pho- 
ton in field m polarized along the angle ips and its entan- 
gled idler twin in field n polarized along the angle tpi at 
the output from the photonic-band-gap structure. 

We are interested only in the part 1^)^ / of the output 

state |V')s"^s,i,¥'i given in Eq. if^ that describes the 
generated photon-pairs. Including the time evolution of 
the free-fields outside the photonic-band-gap structure, 
we can write: 

m))^I = Hffit)) + I'Aff w> + Kf w> + iv'ff w>, 

(27) 
where 
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X a] {iOs)al {uji)exp{iuJst)exp{iuJit)\vac) 



m, n 



F,B- 



(28) 



The creation operators a, 



ot 



SF,iFiSB,iB, intro- 



duced in Eq. (|28|l describe the linearly polarized photons 
outside the photonic-band-gap structure with polariza- 
tion angles (ps and ipi^ and are given as: 



«LK) = Cil)^(^.), C(^^ 






a\^{u,s)^afllS^s), C(c..) = «i"l(c..)- (29) 

We note that we do not explicitly express the dependence 
of quantities on the polarization angles (ps and ipi below. 



C. Properties of the signal and idler fields outside 
the structure 

The mean number iV™" of photon pairs that have a 
signal photon in the frequency interval Aw^ around fre- 
quency ujs and its twin idler photon in frequency interval 
Awi around frequency LOi in field mn is given by 

7v™"(c.„c.,) = (v.rr(t)i^s,„(c.,)n,„(c.,)iV'rr(t)>Ao..Ac.„ 

(30) 
where the operator nj^^{ujj), the density of photons, is 
defined as 

%„> i^3 ) = «]„ i^3 ) djrn. (^J ) ' 3 ^ S,i. (31) 



Using Eq. (PS|l . the expression for iVJ^" in Eq. l(^ can 
be written as: 



T\rmn / \ 



\uJs,UJi)\'^ALUsALU,. (32) 



We get the following expressions for, e.g., the mean 
number iV™"(cJs) of signal photons in frequency interval 
AlOs around frequency Wg'- 

Nritos) = (^™r(i)l".™(^.)I^.T(i))Ac., 

da;,|0™"(a;„u;,)l'AL^,. (33) 



The overall number TV™" of photon pairs emitted into 
field mn is determined by the expression 

/>00 />CXD 

7V™"=/ dios d^,|0™"(^„c^,)|2. (34) 

Jo Jo 

The signal- field energy spectrum S""^{u!s) of field mn 
can be easily determined using the expression 



sric^s) = hio, 



iVI""(c^a) 

Awc 



= TiLJs / rfc^,|(/.™"(w„c^,)l'- (35) 
Jo 

The energy spectra Sj^iujj) {j — s,i; m = F,B) char- 
acterizing the outgoing down-converted fields without an 
inclusion of pair entanglement can be evaluated using 
the signal- field spectra 5™" determined in Eq. (|S5|l : the 
idler-field spectra 5'™" are determined analogously; 



c _ qFF I qFB 

•JSF — JsF + '-'SF ' 

c _ qFF I qBF 

'-Jip — '-'ip ' '-'iF ' 

c _ qBF I qBB 

'-'sB — '-'sB ^ '-'sB ' 

^FB I riBB 



biry J^ 



IB 



(36) 



The properties of the down-converted fields in the time 
domain can be conveniently described using a two-photon 
amplitude Aijs^Ti) giving the probability amplitude of 
detecting a signal photon at time r^ and an idler photon 
at time r,: 



^'""(Ts,T,) 



(vac 



|^W(0,io + T,) 



M+) 



xi?);^(0,to + r.)ICr(M)- (37) 



Assuming the state IV"™") given in Eq. (|28|) . the expres- 
sion in Eq. H37|l for the two-photon amplitude ^™" can 
be rearranged into the form: 



-4""(r„T,) = ^^0™"(r.,T,), (38) 

where the 'Fourier transform' ((/'^'^{ts, Ti) of the function 
(j)^"-"[ujg^uj.i) has been introduced: 



\Ts,Ti} = — / dw, / duj. 



V ^S^l 



X exp{-iujsTs) exp{-iujin). 



(39) 
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FIG. 3: Scheme of Hong-Ou-Mandel interferometer. In a non- 
linear structure a photon from a pump beam p is converted 
into a photon pair copropagating as signal (s) and idler (i) 
beams. After reflection from mirrors M and introduction of 
a mutual delay in a delay line, both photons interfere at a 
beam-splitter BS, and are then detected using a coincidence- 
count detection scheme at detectors Da and Db. 



The photon flux of, e.g., the signal photons, M^"'{ts) 
at time Ts over the transverse profile of area B in field 
mn is defined as 



^rirs) 



eocB 



(cr(^o)i^i:Ho,to 



x^(+)(0,to-l-r,)|V'™r(io)>. (40) 



The photon fiux A/'™"(ts) can be determined using the 
function (^""(ws, 0;^); 
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X 0'""*(w;,w,;)(/>""(ws,t^»)exp(iw;T,)exp(-iw^Ts). (41) 

If the spectrum of the idler field is narrow, we may use 
the alternative expression 



AA" 



\r.) = M 



dr, |0""(T„r,)|^ 



(42) 



Entanglement of the signal and idler photons in a pair 
can be detected in a Hong-Ou-Mandel interferometer (see 
Fig. inj. In order to achieve this type of interference 
between the signal and idler photons, we rotate polar- 
izations of both photons such that they are the same, 
and subsequently introduce a relative time delay ti be- 
tween the two photons. Then two photons arrive at a 
50/50 % beam-splitter whose output ports are monitored 
by detectors. The measured coincidence-count rate Re is 
given by the number of simultaneously detected photons 
at detectors Da and Db, placed at the output ports of 
the beam-splitter in a given time interval. There occurs 
quantum interference between two paths leading to a co- 
incidence count; either a signal photon is detected by 
detector D^ and its twin idler photon by detector D_b or 
vice versa. The normalized coincidence-count rate i?„ in 
this interferometer can be expressed as follows: 



where 

p'""(tO 



2i?J," 



dtA 



dtf 



Re [^'"" {tA ,tB-n )^""* (ts ,tA-Ti)] (44) 



and 



dtA dtB|^'""(iA,ts)r (45) 



-Kn — - 



The symbol Re stands for real part, and the two-photon 
amplitude A™"^ has been introduced in Eq. H37|l . Using 
the expressions in Eqs. (|44|l and (|45|l . we arrive at: 



P""(rO 



Begc2S2 R^ 



-Re 



^0 



(/)™"(ws, cjj)0™"*(wi, cjs) exp(ia;iT;) exp(-Ja;sTi) 



inmn 



(46) 
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(47) 



■K„ (Ti) = l-p (Ti), 



(43) 



D. Cw-limit 

If cw pumping is considered, the spectrum is described 
using Eq. Q). As a consequence the following expressions 
are obtained: 

\<i>-^'-{LO,,U,)\^ = f{u;,,UJ,)5\L0p-L0s-L0,). (48) 

In this case the second power of the Dirac delta function 5 
is formal, so that the above expression should be replaced 
by: 

2T 

|(/)'""(ws,t^»)P = lim --f{L0s,i^i)5{ijJp-u0s-uJ,), (49) 

where the period of nonlinear interaction goes from —T 
to T . Expressions for the physical quantities determined 
above have to be normalized by 2T, which indicates that 
these quantities are related to the period of Is. 



III. TYPICAL CHARACTERISTICS OF THE 

DOWN-CONVERTED FIELDS GENERATED 

FROM A PHOTONIC-BAND-GAP STRUCTURE 

As an example, we study the properties of a nonlinear 
photonic-band-gap structure composed of 25 nonlinear 
layers of GaN of thickness of 117 nm among which there 
are 24 linear layers of AIN with a thickness of 180 nm. 
Material characteristics of GaN and AIN can be found, 
e.g., in p3lll32||. This structure is resonant for the pump 
field at wavelength of 664.5 nm, designed to correspond 



to the first resonance peak near the band edge. This sit- 
uation favors the efficient generation of photon pairs at 
roughly double the pump wavelength, at angles that cor- 
respond to the first, second, and third resonance peaks, 
respectively, for the down-converted fields. We assume 
that the GaN crystallographic x axis coincides with the 
z axis of propagation (see Fig. ^ . In this configuration 
there is strong nonlinear interaction between TE compo- 
nents of the pump, signal, and idler beams. The non- 
linear coefficient d for this geometry is chosen to be 10 
pmV^^. 



A. Cw pumping 
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We first consider pumping by a cw laser tuned at first 
resonance peak, i.e. Xp = 664.5 nm. A photon pair 
whose signal photon is emitted along the angle 9s may 
be described by a two-photon amplitude A with typical 
shapes both in time and spectral domains. In the spectral 
domain, the two-photon amplitude ^(wg, Wi) can be writ- 
ten in the form s{u!s)S{ujp — ujg — Wi), where the function 
s{ujs) is linearly proportional to the spectral amplitude 
of the signal field. This is a consequence of stationarity 
of the emitted down-converted fields. The probability 



b) 



\Airs 



of detecting a signal photon at time Ts and 



an idler photon at time t^ is shown in Fig. ^ for mode 
FF (both photons exit the structure at z — zn). 

The shape of the probability \A{Ts,Ti)\'^ can be easily 
understood when we consider the fact that both photons 
in a pair are born together at the same instant of time, 
and then propagate independently inside the structure. 
This independent zig-zag propagation increases the aver- 
age time delay between the two photons, which takes into 
account all possible realizations of the random process of 
propagation of a photon that undergoes multiple bounces 
inside the structure. The 'random' zig-zag propagation 
causes interference between different paths. This leads to 
a typical structure with local minima and maxima along 
the direction Ts + u ~ const; as shown in Fig.^Ja. Also, 
the greater the difference Tg — u the smaller the value 
of the two-photon amplitude A because the propagation 
with more zig-zags is less probable. We note that typical 
stripes along the direction Ts — Ti — const reflect station- 
arity of the process. 

Typical spectral properties of this structure are docu- 
mented in Fig. |5l where the signal-field energy spectrum 
Sg^ is depicted as a function of the angle 9s for field 
FF. The spectrum Sf^ giving the probability of emis- 
sion of a signal photon at frequency ojs along the angle 
9s is a complex function in its variables because it is 
built up by interference of fields coming from the multi- 
layer structure. A maximum photon-pair generation rate 
is obtained for 6s ~ 13.8 deg for degenerate frequencies 
(ujs ~ tOi). Photons from a generated photon pair can 
exit the structure also at z — zq] their energy spectra 
Sf^{ujs,9s), Sf''{ujs,es) and Sf^iujs,9s) are similar to 
that shown in Fig. |S1 This is because efficient genera- 
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FIG. 4: Probability |^^^p of detecting a signal photon at 
time Ts and its twin idler plioton at time Ti (a) in mode FF 
for cw pumping is sliown; cut of the graph along the line 
Ts + Ti — 450 fs is shown in (b); the probability |^^^p 
is normalized such that one photon pair is emitted, i.e. 
J dTs J dTi\A^^{Ts,Ti)\^ — 1: 9s — 14 deg; logarithmic scale 
on z axis is used. 



tion of photon pairs occurs provided both the signal and 
idler fields are resonant in the structure so that the gen- 
erated photons exit the structure at z = zq or z = zjv 
with comparable probabilities. The signal-field intensity 
transmission |Tsp as a function of normalized frequency 
2ujs/ujp and angle 9s is shown in Fig. for comparison. 

The angle-dependence of the signal-field FWHM (full- 
width at half maximum) AA^ of the energy spectrum in 
wavelengths and value of energy spectrum S"'^'^^ at the 
central frequency ojg are shown in Fig. |7| High photon- 
pair generation rates are observed for the frequency- 
degenerate case {uj'^ fa wf ) for values of the angle 6s ap- 
proximately equal to 14, 35, and 55 degrees. For these 
values of the angle 6s the peaks of the signal-field in- 
tensity transmission |Tsp cross the value ujs equal to 
ujp/2 (see Fig. EJ. The first and highest peaks in the 
signal-field energy spectra Ss (see Fig. [TId) correspond 
to the case when both the signal and idler photons are 
tuned at the first resonance near the band edge. At 
this frequency, localization of the down-converted fields 
is maximum, and the nonlinear process is strongest. The 
other two sets of peaks in the signal-field energy spec- 
tra Ss correspond to the second and third peaks in the 
intensity transmission \Ts\'^. At these frequencies, the 
nonlinear interaction is somewhat weaker due to smaller 
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FIG. 5: Energy spectrum S^^ of mode FF versus normalized 
signal- field frequency 2(jJs/(jJp and angle 9s of emission of a 
signal photon is shown in (a); cuts along the lines 9s — 14 deg 
(solid curve) and 9s = 28 deg (solid curve with *) are shown 
in (b); cw pumping is assumed; logarithmic scale is used on 
the z axis. 
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FIG. 7: FWHM A\s (a) and maximum value S™'"' at the 
central frequency cjj (b) for the signal-field energy spectrum 
Ss{lus) of modes FF (*), FB (A), BF (o), and BB (o) as 
functions of the angle Os of emission of a signal photon are 
shown; cw pumping is assumed. 
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FIG. 6: Signal-field intensity transmission \Fs^ as a function 
of normalized frequency 2lJs/i^p and angle 9s- Positions of 
peaks of resonance as they depend on the angle 9s are clearly 
visible. 



field intensities, i.e., slightly worse localization of the 
down-converted fields. The spectral FWHMs for these 
frequency-degenerate emissions lye in the interval from 
10 nm to 15 nm (see Fig. [7^). For these values of the 
angle Qg the best available constructive interference in 
the structure occurs. In a relatively narrow frequency 
interval the signal and idler fields are enhanced by con- 
structive interference. Outside this frequency interval. 



destructive interference occurs. For the other values of 
the angle ds constructive interference is less pronounced, 
the range of frequencies having constructive interference 
is broader, and spectral shapes with several peaks may 
occur (see Fig-ISJa). We thus obtain down-converted fields 
with spectral FWHMs reaching even 80nm (see Fig. [7^) . 
In this case the FWHM indicates the entire active width, 
which might include several peaks. However, photon-pair 
generation rates are low in this case. 

Quantum correlations (entanglement) between the sig- 
nal and idler photons in a photon pair are visible in the 
profile of the normalized coincidence-count rate -R„(r/) in 
a Hong-Ou-Mandel interferometer (see Fig. [S] for mode 
FF^ . For values of the angle 0^ with strong constructive 
interference [Bg k, 14, 35, 55 deg) this profile is given by 
a relatively broad dip (see Fig.|SlD). On the other hand, 
oscillating curves with several local minima and maxima 
characterize profiles of the normalized coincidence-count 
rates -R„(r/) for other values of the angle Bs- A typical 
profile for this case is shown in Fig.jSla. These profiles re- 
flect broader spectra of the down-converted fields along 
these angles as well as a nonzero difference of the cen- 
tral frequencies w^ and wf of the signal and idler fields. 
Period of these oscillations is proportional to l/(a-'s — tjf). 

Position rf of the central dip, its FWHM At;, and vis- 
ibility V of the normalized coincidence-count rate i?„ as 
functions of the angle B^ are shown in Fig. |51 We can 
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FIG. 8: (a) Normalized coincidence-count rate R^^ of mode 
FF in Hong-Ou-Mandel interferometer as a function of rel- 
ative time delay ti and angle 6s of emission of a signal pho- 
ton; (b) cuts along the lines 6s = 14 deg (solid curve) and 
9s — 28 deg (solid curve with *); cw pumping is assumed. 



see in Fig. |2K that the maximum overlap of the signal 
and idler photon fields is reached for a nonzero mutual 
time delay t; for nonsymmetric modes FB and BF. This 
means that the signal and idler photons exit the struc- 
ture with a small mutual time delay. This nonzero mu- 
tual time delay means that interference effects inside the 
structure prefer the generation of photon pairs from non- 
linear layers positioned closer to one edge of the struc- 
ture, or in any case approximately co-located. The vis- 
ibility V of the coincidence-count rate i?„ equal to one 
is naturally observed for symmetric modes FF and BB. 
Values of the visibility V better than 0.9 can also be 
reached in nonsymmetric modes FB and BF provided 
that the fields are generated along the angles 9s in the 
vicinity of those with a strong constructive interference. 
The structure of quantum correlations between the sig- 
nal and idler fields in time domain is complex and pre- 
vents the visibility V from reaching higher values along 
the angles Og, where there are broad spectra of the down- 
converted fields. The width At; of the dip in the normal- 
ized coincidence-count rate Rniv) determines the entan- 
glement time, i.e. the time interval in which both pho- 
tons can be detected. The FWHM Ar/ of the central 
dip for the structure ranges from 50 fs for values of the 
angle 9s with broad spectra of the down-converted fields, 
to 150-250 fs for values of the angle 9s at where strong 
constructive interference occurs (see Fig.|3;). 
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FIG. 9: Position rf of the central dip (a), its visibility V 
[V = j^] (b), and FWHM Ar; of the dip (c) for the nor- 
malized coincidence-count rate R„{ti) in Hong-Ou-Mandel in- 
terferometer in modes FF (*), FB (A), BF (o), and BB (o) 
as functions of the angle 6s of emission of a signal photon; cw 
pumping is assumed. 



B. Pumping by an ultrashort pulse 

We now assume that the structure is pumped by an 
ultrashort pulse with an unchirped Gaussian profile [see 
Eq.ll^J] having pulse duration Tp equal to 200 fs, and cen- 
tral (carrier) frequency u>p at the wavelength of 664.5 nm. 
The squared modulus |£^pp of the incident pump- field 
amplitude spectrum 'Epp^Tsi^p) fits well into the peak 
of resonance of the pump-field intensity transmission 
\Tp{LUp)\'^, as depicted in Fig. ^ 

Only signal and idler photons with frequencies uJs and 
LUi for which the sum ujs + lu^ lyes within the pump- 
pulse spectrum may be generated. The allowed values 
of the frequencies lUs and uji are indicated in the graph 
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FIG. 10: Pump-field intensity transmission \Tp'^ (solid curve) 
and squared modulus |Epjj,tb| of the pump-field amplitude 
spectrum (solid curve with *) as functions of the normalized 
pump-field frequency uj^/uj'^. 
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FIG. 11: The probability |(;/)^^|^ of having a photon pair 
in mode FF with a signal photon at the normalized fre- 
quency 2uia/ui'^ and an idler photon at the normalized fre- 
quency 2aji/ajp for pulsed pumping; the probability 10^^ ^ 
is normalized such that one photon pair is emitted, i.e. 
A ^ diUs ^ du,\<j)''^ {uj,,LOi)\'{Ljlf = 1; e, = 14 deg. 



of the probability \(j){ujs,ijji)\'^ of emitting a signal pho- 
ton at frequency tOg and its twin at frequency Ui. A 
typical shape of the probability \(j)(ujs,i^i)\'^ is shown in 
Fig. ^ for mode FF. We note that the shape of the 
probability |(/)(ws,Wi)p in Fig. 1111 may be considered as 
being composed of many curves defined above the lines 
LOg + uji — Ljp — const, and it may be assumed to be lin- 
early proportional to those characterizing cw pumping at 
the frequency Up. The probability |^(r5,ri)p of detect- 
ing a signal photon at time r^ and its twin at time t^ de- 
picted in Fig.ll2lshows that both signal and idler photons 
occur in sharp time windows as a consequence of pump- 
ing by an ultrashort pulse. The probability |yl(Ts,Ti)p 
has a typical droplet shape that reflects the fact that the 
longer the photons from a pair propagate inside the struc- 
ture, the greater the average difference of the occurrence 
times Ts and ti of the photons. The structure with local 
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FIG. 12: The probability \A''^\'^ of detecting a signal 
photon at time r^ and its twin idler photon at time Ti 
in mode FF for pulsed pumping; the probability |^^^p 
is normalized such that one photon pair is emitted, i.e. 
i dTs j dTi\A^^ {Ts,Ti)'^ — \\ 6s — 14 deg; logarithmic scale 
is used on the z axis. 



minima and maxima farther from the diagonal in Fig. 1121 
reflects interference in the layered structure. 

Because the pump-pulse amplitude spectrum Ep^^TB 
fits well into the peak of resonance of the structure, and 
the emitted frequencies of the down-converted fields have 
also to fit into the structure, the energy spectra of the 
signal and idler fields with pulsed pumping are very sim- 
ilar to those obtained for cw pumping. Also, the be- 
havior of photon pairs generated by pulsed pumping in 
the Hong-Ou-Mandel interferometer is similar to that ob- 
served with cw pumping; only the visibilities V of the 
normalized coincidence-count rates Rn{Ti) are slightly 
worse. 

The signal and idler fields are now emitted in the form 
of ultrashort, pulsed fields in multimode Fock states. A 
typical dependence of the signal-field photon flux Nf^ 
at time Tg and angle 9s is shown in Fig. ^J Fig. ^1 sug- 
gests that the photon flux Ng typically spreads to longer 
times as a consequence of zig-zag movement, i.e. mul- 
tiple scattering, of photons inside the structure. Time 
delay t^ and FWHM At^ of the pulsed signal held as 
well as photon flux iV™^'^ at the center of the pulsed fleld 
as they depend on the angle 0g are depicted in Fig. E| 
The time delay r^ of the signal pulse with respect to the 
pump pulse ranges from 70 fs to 200 fs. It reaches its 
maximum values for values of the angle 9g with strong 
constructive interference {9g « 14,35,55 deg). FWHM 
Ats of the signal pulse goes from 270 fs to 320 fs. The 
narrower the spectrum of the signal fleld (given by AAg 
in wavelengths) the greater the signal-field pulse dura- 
tion Ats (compare Figs. [7^ and [Hb) . The largest value 
of the photon flux iV™^'^ can be reached for the angle 9g 
equal to 13.8 deg (see Fig. ITlh). 

These results clearly show that the considered 
photonic-band-gap structure maintains a pulsed char- 
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FIG. 13: Photon flux AT^^ of mode FF as a function of time 
Ts of detection of a signal photon and angle 6s of emission of a 
signal photon; pulsed pumping is assumed; logarithmic scale 
is used on the z axis. 



acter of the nonlinear process and allows generation of 
pulsed down-converted fields with an ultrashort dura- 
tion. Such pulsed fields are required in many experi- 
ments with photon pairs in which time synchronization 
of several photon pairs is necessary. 



C. Efficiency of the nonlinear process 

The model that we have developed is one dimensional 
with respect to the spatial coordinates, and does not in- 
clude all aspects (especially those related to transverse 
profiles of the interacting optical fields) of the nonlinear 
process in a real structure. For this reason we do not 
determine the absolute values of photon-pair generation 
rates. Instead, we judge the efficiency of the suggested 
structure (stemming from constructive interference of the 
interacting optical fields) with respect to an ideal refer- 
ence structure which fully exploits the nonlinearity, but 
does not rely on interference. 

This reference structure has formally all linear indices 
of refraction equal to one, so effects on the boundaries 
between layers are suppressed. It is also assumed that 
the nonlinear process is fully phase matched. The ori- 
entations of nonlinear layers, as well as polarizations of 
the interacting optical fields, are such that the greatest 
nonlinear effect occurs. A photon pair emitted from this 
structure is described by the following output state 1^)1^1 
[compare Eq. (^0)) ]: 
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(50) 
where a^(a;s) [a]^(a;i)] denotes a creation operator of a 
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FIG. 14: Time delay r^ (a), FWHM Ar, (b), and photon 
flux A*'™^" at the center (r^ — rf ) for the pulsed signal field in 
modes FF (*), FB (A), BF (o), and BB (o) as they depend 
on the angle 6s of emission of a signal photon; pulsed pumping 
is assumed. 



photon in the signal [idler] field. The function max used 
in Eq. (|50|l gives the maximum value from the elements 
of a tensor in its argument. 

The relative photon-pair generation rate 77™" (ws) of a 
pair with a signal photon at frequency cOs in mode mn is 
then given as: 



^™(^.) 






(51) 



where the signal-field energy spectrum S*™" is given in 
Eq. PB|| , and the signal- field energy spectrum Sl^^ of the 
reference structure is determined in the same way using 
the output state \ipY°l written in Eq. (|50|l . 

The relative photon-pair generation rates rj^^, V^^^ 

-°/2 



^BF 



and 77 for the signal-field frequency uig 
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FIG. 15: Relative photon-pair generation rates rj {*), V 
(A), rj^^ (o), ri^^ (o), and overall relative photon-pair gen- 
eration rate rj {rj = f)^^ + rj^^ +rj'^^ +ri^^) at the signal-field 
frequency Ws = ujp/2 for cw pumping as a function of angle 
^3 of emission of a signal photon. 



as a function of the angle Os are shown in Fig. ^] as- 
suming cw pumping. The relative photon-pair generation 
rates for pulsed pumping are lower in comparison with cw 
pumping, because frequencies at the edges of the pump- 
field spectrum have lower efficiencies than those in the 
middle. Despite this, the considered pulsed pumping has 
the overall relative photon-pair generation rate rj equal 
to 9, whereas ry = 12.5 for cw pumping for 6's « 13.8 deg. 
The structure that we consider is approximately 12 times 
more efficient than the reference structure, provided we 
consider all the emitted photon pairs (given by the overall 
relative photon-pair generation rate 77 in Fig. ll5() . Even 
when we restrict ourselves only to the mode FF^ we still 
have an enhancement of the nonlinear process by a factor 
of three. For comparison purposes, one layer of GaN of 
thickness 25 x 117 nm (i.e., containing the same amount 
of nonlinear material as our structure) has the overall 
relative photon-pair generation rate 77 equal to 0.09 for 
values of the angle 9s less than 20 deg. The inclusion of 
linear layers of AIN inside the structure and the accom- 
panying interference effects thus increase the photon-pair 



generation rates by two orders of magnitude. 



IV. CONCLUSION 

We have developed a quantum model of spontaneous 
parametric down-conversion (generating entangled pho- 
ton pairs) in one-dimensional, nonlinear, layered media. 
Using the model we have determined a two-photon ampli- 
tude, and we have provided measurable characteristics of 
the down-converted fields: marginal signal and idler en- 
ergy spectra, time-dependent photon fluxes of the signal 
and idler fields, coincidence-count interference patterns 
in the Hong-Ou-Mandcl interferometer, and photon-pair 
generation rates. 

A specific structure was suggested as an efficient source 
of photon pairs, and interference effects of the interact- 
ing fields can enhance the photon-pair generation rates 
by as much as two hundred times. The widths of the 
spectra of the down-converted fields and entanglement 
time of photons in a pair depend strongly on the angle of 
emission, and vary from 10 nm to 80 nm and from 50 fs 
to 250 fs, respectively. Pumping the structure with an 
ultrashort pulse of time duration of several hundreds of 
fs can lead to the generation of pulsed signal and idler 
fields extending over several hundreds of fs. 

Therefore, we conclude that nonlinear one-dimensional 
photonic-band-gap structures represent a promising, new 
efficient source of photon pairs. 
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